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Abstract 



^ ■ A new method for calculation of traces of Dirac 7-matrices contracted with massless 

, vectors in Minkowski space is discussed 

o 

^ 1 Introduction 

^ ■ Calculation of Dirac 7-matrices traces is standard procedure in high energy physics computa- 
tion. However the expressions for the traces of more than ten 7-matrices are too long. 
" Application of Chisholm pj identities and Kahane [2J algorithm leads to short expressions 
^ \ for traces in 4-spacetime. More compact formulae can be obtained through the new method 
' proposed in [3] and based on the properties of orthonormal bases in the Minkowski space and 
isotropic tetrads (see e.g. [l]-[5]) constructed from the vectors of these bases. 

The purpose of the article is further simplification of the formulae presented in [3] for the 
case of Dirac 7-matrices contracted with massless vectors. This task is important today because 
most of analytical calculations in high energy physics are performed at massless approximation 
(see, for example, [6]). 

The article is arranged in the following way. The new method of trace computation for 
arbitrary vectors is shortly described in Section 2. In Section 3 one explores application of the 
traces calculation method for Dirac 7-matrices contracted with massless vectors. The formulae 
obtained there are good enough to realize an algorithm for any numerical calculation. 

Note, that the method of calculation of the traces from the [3] is included in the system of 
analytical computation ALHEP [8] , and this article contains the proposals for further improve- 
ment of the similar programs. 
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In this paper we use the Feynman metrics: 

/X = 0, 1, 2, 3, a'^ = (ao, a), a^ = (oq, —a), {ab) = a^b^ = ao^o ~ (^b , 
and a sign of the Levi-Civita tensor is fixed as 

^0123 = +1 • 

Orientation of the orthonormal basis vectors in Minkowski space 

{IaIb)=9ab, (A 5 = 0,1, 2, 3) 

is constrained by 

i/J-ivipiX _ _|_-| 
t^i/pA''0 ''1 ''2''3 — "T-L . 

Isotropic tetrads q±, e± are constructing in the following way: 

e+ = , e_ = ^ , (e+e_) = -1 , = e± ; ^ ^ 

(g+e+) = (g+e_) = (g_e+) = (g_e_) = . 

2 The new method of trace calculation for Dirac 7-matrices 
contracted with arbitrary vectors 

The new method to calculate traces of Dirac 7-matrices contracted with arbitrary vectors 
a, (i = l,---,2n) in Minkowski space is presented in [3j. So the following expressions were 
obtained: ^ 

- Tr[(l - 7^)0102] = Fi(ai, as) + ^3(^1, «2) ; (2) 

^Tr[(l - 7^)01020304] = 

= Fi(oi, 02)^1(03, 04) + ^2(01, 02)^4(03, 04) + ^3(01, 02)^3(03, 04) + ^4(01, 02)^2(03, 04) ; 

(3) 

where Fi, F2, F3, F4 are some functions (see below) of vectors Oj. 
Hawing calculated 

Tr[(l - 7^)01 ■ ■ ■02„_ia2„] , 

we may obtain an expression for 

Tr[(l - 7^)01 ■ ■ ■ 02„_ia2„a2„+ia2„+2] 



from the previous one by the following replacement: 

Fl{a2n-1, 0,271) Fi(a2n-1, (''2n)Fi{a2n+l, a2n+2) + -^2(a2n-lj fl2n)-f4(«2n+l, 0,2n+2) 

F3{a2n-1, 02n) -^3(o^2n-l) Q2n)-^3(Q2n+l, 020+2) + -^4(0271-!, 0^2n)-^2(o^2n+l) Q2n+2) 

F2{a2n-1, 0,2n) Fi[a2n-1, Ot2n)-p2(ot2n+l, 0t2n+2) + -^2(a2n-lj ti2n)-^3(ti2n+l, 0t2n+2) 

F4{a2n-1, 0,2n) F^{a2n-1, Q2n)-^4(o2n+l > Q2n+2) + -^^4(0271-1, 0'2n) Fi{a2n+1 , Q2n+2) 

The expression for 

Tr[(l - 7^)61 • • • a2n-ia2n\ 
obtained through the method will contain 2" terms. 

Tr[(l + 7^)aia2 • • • a2n] = (Tr[(l - 7^)0102 • • • a2n]) * • 
The functions mentioned above have the following forms: 

Fi{ai,aj) = 2[{aiq-){ajq+) - (0^6+) (%€_)] = 



1 1 
= 1 Tr[(l - j^)q+q-aiaj] = -- Tr[(l - j^)e_e+aidj] , 

F^{ai,aj) = 2[{aiq+){ajq-) - {aieJ){aje+)\ = 

= \ Tr[(l - l^)q-q+aiaj\ = Tr[(l - j^)e+e-aidj] , 



F2{ai,aj) = 2[{aie+){ajq-) - {aiq-){aje+)] = 2 G ^ 



e+ g_ 



= -Gr^ ^0 +^Gr^ '^n + G -^^r^ 

V to t2 / V ti t3 / \n h J V to t3 

= ^1^[(l-7')g-e+aia,-] , 



F4(ai,aj) = 2[(aie_)(ojg+) - (0^5+) (a^e-)] = 2 G 



= - G ( a, N ^ ^ / ai aj\ _ ^/ai aj \ ^ ■ ^ f ai % 
V to t2 / V ti t3 / Kk h J V to t3 

= ^ Tr[(l - 75)g+e_Oiaj] , 



where G are Gram determinants. 

There is an example of Fj functions for some orthonormal basis. Let's to fix the last as 
follow: 



/o^ = (l, 0,0,0), /r = (0, 1,0,0), i2^ = (0, 0,1,0), /^ = (0, 0,0,1) 
Then the isotropic tetrads can be expressed as 

g± = ;^(1, ±1,0,0), (ag±) = -^(aoTOa;) , 
e± = ;^(0, 0, 1, ±i) , (ae±) = -;^(«2/ ± ^a^) • 

At last 



^l(«i,«i) = 


-- {ttittj) - [ {ai)o{aj)^ - 




] + {ai)y{aj), - 


iai)z{aj)y ] , 


F3(oj,aj) = 


-- {aittj) + [ {ai)o{aj)^ - 


{ai)x{0'j)o 


iai)y{aj), - 


iai)ziaj)y ] , 


F2 


[ai,aj) = [ {ai)o{aj)y- 


{ai)y{aj)o ] 


+ i[ {ai)x{aj)^ - 


{ai)z{a'j)x ] + 




+ [ iai)x{aj)y - 


{ai)y{aj)x 


] + {ai)o{aj)z - 


{O'i)z{0'j)o ] , 




'fluttj) = [ {ai)o{aj)y - 


{ai)y{aj)o ; 


+ i [ {ai)x{aj)z - 


{ai)z{aj)x ] - 




- [ iai)x{aj)y - 


{ai)y{aj)x 


(ai)oK)z- 


{ai)z{aj)Q ] . 



For clarity one can writes expressions of traces in the following form: 



Tr 



(1 - 7^)01620304 • • • a2„_ia2ri 



(10) 



(11) 



(12) 



Tr 



-^3(01,02) ^^4(01,02) 
^2(01, 02) Fi(ai,a2) 



^3(03, 04) ^4(03, 04) 
^2(03, 04) -^1(03,04) 



F3{a2n-1, tt2n) F4^{a2n-1 , tt2n) 
F2{a2n-l,0'2n) Fi(a2n-1, 1271) 

(13) 

It's obvious that the trace of 2n Dirac 7-matrices is reduced to trace of n matrices with 2x2 
dimension. 

For calculation of expressions like this 

Tr (7^0162 • • • a2n+i) Tr(7''6i62 • • • fe2m+i) , 
where summing over index p is supposed, one can use the Fiertz transform: 

[(1 ± 7')7.]i.[(l T l')lp]ki = 2(1 ± T 7^),,- . (14) 
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So one achieves the following expressions: 

Tr[(l ± 7')7paia2 ■ ■ ■ a^n+i] ■ Tr[(l T l^'h'kk ■ ■ ■ Kn+i] 
= 4 Tr[(l =F 7^)aia2 ■ ■ ■ a2n+i&i^2 ■■■ h 



(15) 



Tr[(l ± 7')7paia2 ■ ■ ■ a2„+i] ■ Tr[(l ± 7')7"^i&2 ■ ■ ■ km+i] = 

= Tr[(l ± 7')7pai«2 ■ ■ ■ a2n+i] ■ Tr[(l T l^'h'hm+i ■ ■ ■ kk] = (16) 

= 4 Tr[(l =F 7^)01^2 ■ ■ ■ a2n+i^2m+i ■ ■ ■ hh] , 

Tr(7paia2 ■ ■ ■ a2„,+i) ■ Tr(7''6i62 ■ • ■ 62m+i) = 2 Tr \a1a2 ■ ■ ■ a2n+i (hh ■ ■ ■ hm+i + hm+i ■ ■ ■ hh) 

(17) 

3 Trace calculation of Dirac 7-matrices contracted with 
massless vectors 

3.1 Application of massless vectors 

In the method of trace calculation, what was briefly explained in Section 2, vectors contracted 
with Dirac 7-matrices are arbitrary. But there is a significant simplification of expressions for 
traces when vectors are massless (see e.g. [6], [7]). 
For massless vectors the particular equation is true: 

a? = {a,lof - {aihf - {a,hf - {adzf = , (18) 



I.e. 



or 



{ailof - {aihf = {aikf + {a^hf , (19) 

(aig+)(aig_) = (0^6+) (a^e.) . (20) 

Thus the functions Fk{ai,aj) become to be like the function Fi{ai,aj) multiplied by some 
factor: 

F2{ai,aj) = 2[{aie+){ajq-) - {aiq-){aje+)] = 

= 2[(a,e+)^^^^±l^^ - (a,g_)(a,e+)^] = -|^2[(a,g_)(a,5+) - (a,e+)(a,e_)] , 

(21) 
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that is 

F2{ai,aj) = -j^^^^Fi{ai,aj) . 
[ajq+j 

In the similar way 

F(n n\- F (n n \ - i^^l^ F ( n n\- i^^l^ F ( n n\ 

(aig_) {ajq+) (aig_) {ajq+) 

Fi{ai,aj) = Fi{ai,aj) . 

Formulae of the traces become to be simpler by far. 
Through the identity 

(1 ± 7')gg(l ± 7')g = Tr[(l ± ^')qQ] (1 ± f)q , 

which is valid for any massless vector q and any operator Q, one can obtain 

1 1 

- Tr[(l - 7^)aia2a3 ■ ■ ■ a2„] = — Tr[(l + 7^)^_aia2a3 ■ ■ ■ a2„ai] = 

2 4(aig_) 

1 1 

Tr[(l + 7^ 010203(1 + ■ ■ ■ a2nai] = 



4(aig_) 4(a3g- 
Tr[(l + 7^)^-010203] 1 



Tr[(l +7^)9-%- ■■«2nai] 



4(oig_) 4(o3g_ 
Tr[(l + 7^)^-010203] Tr[(l + 7^)^-030405] Tr[(l + 7^)g_a2„-i02nai] 



4(oig_) 4(a3g_) 4(o2„-ig-) 

Further 

4(^ Tr[(l + 7^)g-0.o,o.] = ^^^^2){a,q^) ^'^^^ + ' ^')^^-M 

Tr[(l + 7^)^-aiOj^+(l + 7^)^-g+ajafc] = 



32(oig_)(ojg+) 



1 rp r/. — irp , . . . i Fi{ai,aj)F^{aj,ak) 

Tr[(l - 7=')g+g-OiOjJ Tr[(l + 7'')g-g+OjOfcJ - 



32(oig_)(ojg+) 2(oig_)(ojg+) 
and finally from (EEl), (ET]) it follows that 

- Tr[(l — 7^)0 O O a 1 _ -^ll^^li '^2)-F3*(Q-2i Q3) • • • -^l('^2n-l; a2re)-^3*(Q2n; Ql) 

2 1 2 3 2n 2"(0ig-)(02g+)(03g_) ■ ■ ■ (02ng+) 

From ([22]), ([23D we have 

Fi(oj,Oj)F3(oj,Ofc) _ F2{ai,aj)F2(aj,ak) 
(ojg+) (oj-g-) 
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and (1251) takes the form 



1 rp r/1 - - 1 _ ^ 1^r^ -^2 («1 , 0-2) -^2 (^2 , 0.3) " " " -^2 (o.2n-l , 02n) ^2* (^2n j «,! ) 

-Tr[(l-7)a,a2a3---a2.]-(-l) 2Ha,q.)ia,q.){a,g.) ■ ■ ■ {a,^g,) " 

(30) 

Note that 

Tr[(l - 7^)010203 ■ ■ • a2n] = Tr[(l + 7^)0203 ■ ■ ■ a2n«i] = (Tr[(l - 7^)0203 ■ ■ ■ 02^01]) , (31) 
then (|28l) leads to the one more expression of traces 



-Tr[(l - 7^)ai«2a3- ■ ■a2„] 



^3(01, 02)^^(02, 03) ■ ■ ■ F3{a2n~i,a2n)F^{a2n,ai) 



2"(aig+)(a2g_)(a3g+) ■ ■ ■ (a2ng- 
At last an identity [see ([23]), 1^] 

F3{ai,aj)Fl{aj,ak) _ Fi{ai,aj)Fl{aj,ak) 



provide for 
1 



Tr[(l-7 )aia2a3- ■ ■a2„ 



^4(01, 02)^4* (a2, 0.3) ■■ ■ Fi{a2n-i-, a2n)Fl{a2n, «i) 



(32) 



(33) 



-1 



2"(ai?+)(«2g+)(a3g+) ■ ■ ■ (a2ng+) 



(34) 

The expressions for traces (128|) . (130!) . (1321) . (l34l) are equivalent. Since in the case of massless 
vectors 

\Fi{ai,aj) \ = 2^J{aiaj){aiq^){ajq+) , 



1-^3(0*, flj) I = 2^{aiaj){aiq+){ajq^) 
\F2{ai,aj) \ = 2^ {aiaj){aiqS){ajqS) 
\Fi{auaj) \ = 2^J {aiaj){aiq+){ajq+) 



(35) 



then 



(see also [6]). 



^Tr[(l - 7^)ai«2a3 ■ ■ ■a2ri 



2"^ (aiOz) (a2a3) ■ ■ ■ (a2nai) , 



(36) 



3.2 Calculation of traces in case of appearance the type ^ uncer- 



tainty 







The appearance of uncertainty type - is possible during numerical calculation by reason of 
denominators presence in the formulae (128!) . (!30|) . (l32l) . 
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For example, when the formula ([31]) is being used and orthonormal basis is chosen in ac- 
cording to ffTUj) . the uncertainty will appear, if 

tti = const ■ q+ 

i.e. if 

The simplest way to avoid the appearing of uncertainties in such points of phase space is 
to use for calculation here another basis vectors, because obtained formulas are correct for any 
orthonormal basis and numerical results received are independent from it's choice. But such 
approach results in unnecessary complicating of the computer program. 

There are three solutions of this problem: 

1. At the points of phase space, where denominators are 0, one should use common formulae 
from Section 2 (see e.g. ^). 

2. One can choose an arbitrary 4- vector t normalized by condition = 1 to perform the 
identical transformation of the initial expression: 



^ Tr[(l - 7^)0-10203 ■ ■ ■a2n] = ^ Tr[(l + '^^)iaiiia2tia^i ■ ■ ■ta2„f| 
= \ Tr[(l + -f^)a[a'^a'^ ■ ■ ■ a'aj = ^ Tr[(l - 7^)^203 " • " ak^i] > 



(37) 



where 
Thus 



a[ = -at + 2{ait)t . (38) 



(O^ = {aiY - 4{aity + Aiaityt' = (a^)^ = , (39) 



so one can use the formulae from Section 3.1 to calculate the transformed expression. At the 
same time the denominators take forms 

(a-fe) = -(««&) + 2(oit)(tg±) , (40) 

where the second term in (l40l) is positive. 

3. It is possible to transform the assumption formula. Let us suppose that is using for 
calculation and for i = 2s 

(o2sg+) = . 
In this situation one can replace [see fl33|) ] 



F4,{a2s-1, Cl2s)Fl{a2s, 0-28+1) -^3(^25-1, 02s)i^i (02^, 02s+ 



(a2sg+j (a2sg- 

Notice that in this case 

{a2sq-) ^ . 



(41) 
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4 Conclusion 

The simple and compact formulae are proposed to calculate traces of Dirac 7-matrices con- 
tracted with massless vectors. These formulae may be easily implemented as a simple yet 
efficient computer algorithm. 
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